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ABSTRACT 

We study polar Alfven oscillations of relativistic stars endowed with a strong global 
poloidal dipolc magnetic field. Here we focus only on the axisymmetric oscillations 
which are studied by evolving numerically the two-dimensional perturbation equations. 
Our study shows that the spectrum of the polar Alfven oscillations is discrete in 
contrast to the spectrum of axial Alfven oscillations which is continuous. We also 
show that the typical fluid modes, such as the / and p modes, are not significantly 
affected by the presence of the strong magnetic field. 
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1 INTRODUCTION 

The soft gamma repeaters (SGRs) are the objects radiating sporadic X-ray and gamma-ray bursts, whose typical luminosities 
are around 10*^ erg/s. Rarely the SGRs were associated with the emission of stronger gamma-rays with peak luminosities 
of about 10'''' ~ 10'*'^ erg/s. These phenomena are typically referred as "giant flares". Up to now at least three giant flares 
have been detected in three different SGRs, which are : the SGR 0526-66 in 1979, the SGR 1900-1-14 in 1998 and the SGR 
1806-20 in 2004. During these giant flares in SGRs we observe an initial strong sharp burst followed by a decaying tail 
which may last for hundreds of seconds. The timing analysis of the decaying tail of the latest two giant flares revealed the 
existence of quasi periodic oscillations (QPOs), whose speciflc frequencies are approximately 18, 26, 30 , 92, 150, 625 and 
1840 Hz for SGR 1806-20, and 28, 53, 84 and 155 Hz for SGR 1900+14 f see I Watts fc Strohmave"il (|2006l ) for a review). The 
most promising mechanism to produce these QPO s are are the magnetar mode ls. Magnetars are thought to be neutron stars 
with strong magnetic fields, £?> 5 x 10^^ Gauss ( Duncan &: Thompson 1992f ). According to this scenario, starquakes may 
be driven during the release of the accumulated magnetic stress in the crust region. The excited seismic oscillations can be 
decomposed into two types, i.e. the axial-type and the polar-type oscillations. Since the axial-type oscillations do not involve 
density variations, it is thought that it will be easier to excite them instead of the polar-type ones. In this case, the magnetic 
field globally and/or the shear modulus in the crust, act as the necessary restoring force for the axial-type oscillations of the 
non-rotating magnetars. 

For the last few years the attempts to e xpla i n the above QPO frequencies in SGRs by using the crust torsional oscillations 
have been based both in N ewtonian (e.g., iLed 1 20071 )1 and in general relativistic dynamics (e.g., ISamuelsson &c Anderss* 
(2007]); ISotani et al.l (|2007l )). These attempts were partially successful, i.e., some of the observed frequencies are in good 



agreement with the frequencies of torsional oscillations. But very soon it has been realized the difficulty to explain all the 
observed frequencies in SGRs by using only the torsional oscillati ons of the cru st (magnetized or not). In an attempt to explain 
the QPOs without making use of the crust torsional oscillations. Levin ( 20061 ) using a toy model suggested that the torsional 
Alfven oscillations could not form a discrete spectrum b ut instead they should for m a continuum. He strongly pointed out 
the importance for the crust-core coupling and soon later iGlampedakis et al. I liooi) came with a similar suggestion i.e. that 
the QPO spectrum may be explained by considering the torsional Alfv en oscil l ations of the magnetar' s core. Abou t a year 
later, this suggestion has been verified with a more elaborate model by Levin (2007), while next year Sotani et al. ( 2008al ) 
working on a more realistic model showed that the spectrum of the torsional Alfven oscillations of relativistic magnetars 
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forms a continuum. Moreover, they demonstrated that there exist two distinct families as the so called "lower" and "upper" 
QPOs. Finally, they figured out that the observed QPO frequencies in SGRs (at least the lower observed frequencies) can be 
explained by using the frequencies of torsional Alfven oscillations as well as the crust torsional oscillations. 



This la st result about the torsiona l Alfven oscillations, has been studied in more detail in two recent papers bv lColaiuda et al 



and ICerda-Duran et al.l I l|2009h . One could actually argue that the features of torsional Alfven oscillations are quite 
well understood. However, the spectrum of the polar Alfven oscillations of magnetars is quite unexplored. So far, there are 
only a few studies of polar oscilla tions of the magnetized stars, where for example, they consider only the fluid modes for Ap 
stars dShibahashi fc Takata 19931) . the oscillations in a spherical shell for Newtonian stellar models wit h an incompressible 
fluid (jRincon fc RieutordlboO^ and the oscillations in the crust region of Newtonian models (|Leg|2007|). 

From the previous discussion one can realize that the polar part of the spectrum of relativistic magnetized stars has not 
yet been studied in detail while it is still unknown whether spectrum of the polar Alfven oscillations is discrete or continuous. 
Moreover, some of the higher QPO frequencies observed may be related to polar and not to axial Alfven oscillations. Finally, 
polar oscillations are related to the emission of gravitational waves and the ultimate result will be the e stimation of the 
amount of the energy released in the form of gravitational waves during the giant flares ([Abbott et al.ll2007r ). 

In an attempt to shed some light to the various features of polar Alfven oscillations, we performed two-dimensional 
numerical evolutions of the perturbation equation of relativistic magnetar models endowed with a global poloidal dipole 
magnetic field. It should be mentioned here that in the present study we omit the effects due to the presence of a crust and we 
focus only on the axisymmetric polar oscillations of the core. The omission of the crust simplifies considerably the boundary 
conditions of the problem without affecting significantly the actual frequencies of the oscillation modes. As a future project we 



can fo resee the inclusion of the crust which induces additional families of shear and interfacial modes (e.g.. IVavoulidis et al 



1 2008t )) which can be used to explain the observed QPO frequencies. In other words, this article is the first step towards a 
detailed study of the oscillation spectrum of magnetized relativistic stars. 

This paper is structured as follows. In the next section we describe the equilibrium configuration adopted in this paper. 
In section [3] we drive the basic perturbation equations with the appropriate boundary conditions. The numerical results are 
shown in section 14] and finally we give a conclusion in section (5] Unless otherwise noted, we adopt units of c = G = 1, where 
c and G denote the speed of light and the gravitational constant, respectively, while the metric signature is (— , +, +, +). 



2 EQUILIBRIUM CONFIGURATION 

We assume that the background models of the strongly magnetized stars that we will study, are spherically symmetric and 
non-rotating. This choice is justified since all known magnetars rotate very slowly, i.e., periods of a few seconds while it 
is true that even for magnetars with stro ng magnetic fields suc h as B > 10^^ Gauss , the deformation due to the magnetic 
pressure is not significant (see for example IColaiuda et"ZI (|2008l ) ; iHaskell et all (|2008l l ) 0. This is expected since even for the 



case of magnetars the magnetic energy is much smaller than the gravitational binding energy, i.e., £m/£g ~ 10^*(-B/10^^ 
Gauss)^, where £m and £g are the magnetic and gravitational binding energies, respectively. In other words, we consider a 
static spherically symmetric neutron star, which is a solution of the Tolman-Oppenheimer-Volkov (TOV) equations. In this 
case the metric is defined as 



(1) 



where e~^"^ = 1 — 2m{r)/r and the fluid four velocity is — (e~*, 0,0,0). On this stellar model we superimpose a dipole 
magnetic field. We assume that the star consists of a perfect fluid and in addition we make the assumption of the "the 
ideal MHD approximation". This implies that a comoving observer cannot feel the presence of an electric fleld. Then the 
stress-energy tensor is given by 

r''" = (e + p)u''u'' + pg''" + H"Hc, (u''u" + ip^"') - H^H^ (2) 

where e, p, and H'^ are the energy density, the pressure, and the normalized magnetic field defined as H'^ = _B'^/\/47r, 
respectively. In particular, in this paper we adopt the pure poloidal dipole magnetic field, while the effect due to the toroidal 
magnetic fields will be left for future studies. The pure poloidal magnetic field is given as 

{Hr,He,H^)=\ — — ^cosS, — ' sm6',0 , (3) 



^ In this paper, for simplicity, we will not take into account the deformation of the star due to the presence of the strong magnetic field. 
These type of deformations will potentially affect, quantitatively, the oscillation spectrum, although one should not expect any important 
qualitative effect. 
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where ai(r) is a radial component of the electromagnetic four-potential, which is the same as in ISotani et al. (2007). This 
potential ai(r) is actually determined as the solution of the differential equation 



+ — ($ 

dr^ dr 



. , da-i 2 2A , 2A • / ^ 

A)— -e ai = -47re ji(r), 

dr 



(4) 



where ji{r) is a radial component of the four-current given by ji{r) = for^{£+p) where fo is a constant. In order to determine 
the form of ai(r), we impose the regularity condition at the stellar center, i.e., ai{r) — olcT^ + 0(r*), while at the stellar 
surface ai(r) should be smoothly connected to the exterior vacuum solution. 



3 PERTURBATION EQUATIONS 
3.1 Axisymmetric Polar Perturbations 

We restrict attention to axisymmetric polar-type perturbations with the relativistic Cowling approximation, that is, the metric 
perturbations are ignored i.e. we set ^g^i/ = 0. Note that the polar type of perturbations are indepe ndent of the axial-type 
ones, which can be understood by recalling the nature of the axisymmetric perturbations (|Led 120071 ). Then the Lagrangian 
displacement vector for axisymmetric polar-type perturbations can be written as 

e 



= —w[t,r,e), 
= ^\v{t,r,e), 



= 0. 



(5) 

(6) 
(7) 



Ifere we should note that the angular dependence of Lagrangian displacement vector ^* for polar parity can be in principle 
described via a spherical harmonic decomposition, i.e., ^"^ ~ Vim, ~ dgYim, and ~ d,j,Yim/ s\v? 9. Then, the components 
of the perturbed 4-velocity Su'^ in the Cowling approximation can be written as 

= Su*" = 0, 



5u^ = — e 



(8) 
(9) 

(10) 



where dt denotes the partial derivative with respect to t. Additionally, the Lagrangian variation of the baryon number density. 
An, is defined as An/n = — V^'^'^*, where V'^' and A denote respectively the covariant derivative in a 3-dimensional spatial 
part of the metric and the Lagrangian variation. With the earlier definition for the Lagrangian displacement vector, the 
Lagrangian variation of the baryon number density takes the form 

e-^ dW 



An 
n 



dr 



+ 



dV_ 

09 



+ cot 9V 



(11) 



Furthermore, for adiabatic perturbations, with the help of the first law of thermodynamics, the perturbed density and pressure 
are given as 



Sp 



, . An 9e ^„ 



IV- 



An 



n dr ^ 
dp r 



n dr 

where 7 is the adiabatic constant defined as 
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/ainp\ _/Ap\/An\-^ 
\dlnnj s \ p J \ n J 



(12) 
(13) 

(14) 
(15) 



and we also use the following relation between the Lagrangian perturbation A and Eulerian perturbaiton S; 
Af{t,r)c^Sfit,r) + ^C- 

The above assumptions lead to the following for m of the r- and ^-c omponents of the linearized equation of motion with 
vanishing shear modulus, i.e., the equation (35) in lSotani et al.l (|2007l ) with /i = 0, 

2^ 



- (5e 4- 5p) $'e 
+ [-2$'e''''i?^ + H' 



e-2^5p,.+ A' 



2re-^^H'' + cot9H'' 
- 2re-^^ sin" 9H'^ -e 



-2A rr 

e Ho 



(16) 
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$' + A' + - ) ff" + i?",, - -^Hr, 



+ cot OH" 



SH" + H''5H\ + H'^SH" 



,4' 



2 sin 6* cos OH 



1 



SH"' ^ -H^SH%, 



(17) 



where a prime (') denotes the derivative with respect to r, i.e., ' = dr- On the other hand, the Unearized induction equation 
(37) in lSotani et al.l (|2007l ) yields the following relations for the first time-derivative of the components of SH^ 



+ - ) ff" + 



;,| Su'' + (cot eW + H\e) 5u^ + W&u^g - H^vT 
I [<t>' + M ^'^^H" + Su'' + (cot OH'^ - <l>'iy'' + H\g) Su" + ff''^^, - H'^Su^ 

I (^$' + A' + i?^ + i?"^,! ^it" + (cot ei?"* + H'^g) Su" + i?"^ {5u\, + <5?i%) . 



(18) 
(19) 

(20) 

(21) 
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Tv _ 



Substituting the previous relations into the linearized equations of motion p6p and p7[) . one can obtain the following coupled 
system of evolution equations for the two functions describing the fluid perturbations 

(22) 



(23) 
(24) 
(25) 
(26) 



where the coefficients .4oo, .4oi, .4io, and are 
Ao = {e + p + H'Hg + H'f'H^) ^e-^'^-'', 



Aoi = -^e ^"^H^Hg, 

A 1 -2#-A rj- 

J\\Q — T-e ti tiri 



Ax 



(6 + p + R-'Hr + H'^Hi) —e 



while the two expressions Tw and Tv are 



Tw = — (i5e + 5p) $'e '^^ — e ^^5p 



+ [-2$'e-'*i/9 + ff^e - 2Te-^^H'' + cotei/-- - e~^^Hg^r] SH" - e~'"'HgSH\, 
+ [-2$'e"'*i?<^ + ff';^ - 2re-''^ sin^ OH'^' - e-'^^H^,r] SH'>' - e'^^ H^SH'^^r, 



(27) 



Tv 



+ A' 



J/'' + cot OH" 



- Hr 



(Hg 



4> — 



2 sine cos ei?" 



(28) 



Notice that the above system of evolution equations has been written for a general equilibrium magnetic field H^{r, 0, (j)), i.e. 
the form of the magnetic field has not yet been specified. 



3.2 Dipole magnetic field 



As we mentioned earlier, in this article we study the perturbations of dipole fields. In order to simplify the boundary condition 
at the stellar surface, we introduce two new functions w{t,r,6) and v{t,r,9) defined as 



w 

V 



eWsinO, 
eVsinO. 



(29) 
(30) 



Notice that in the new functions the trigonometric term, sin 6, has been introduced in an attempt to improve the numerical 
stability of the code. Then, by restricting our attention only to axisymmetric dipole poloidal magnetic fields describe by (|3]), 
we get a simplified form of the linearized equation of motion (|22l) . that is 



.4oo 

^10 



Aoi 
All 









Tw 






= esin^ 










Tv _ 



(31) 
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where the coefficients .4oo, .4oi, .4io, and become 

-2S-A 



47rr 



1 / _A / . 

e ai siny) 



. aiai _2'S-2A 

.4oi = — -e sin y cost 



^10 = 

^11 = 



27rr-2 
aifli' -2#-A 

-2* 



sm f cos f 



(e + p)r H cos ( 

7rr^ 



and the expressions J-'w and will be simplified to the following forms 



J^w = — {5e + 5p) $ e — e 5p,r + 



J^v = —Sp.e + V4^e^jiSH^ sinO + - — 



(32) 
(33) 
(34) 
(35) 

sin 61, (36) 
(37) 



Notice that in order to derive the above equations, we used the relation 
A' + -) SH'' + 5H\,. + 5H'^,e + cot OSH" = 0. 



(38) 

This relation is satisfied throughout the star for the axisymmetric polar perturbations. 

The components of the perturbed magnetic field 5H'' defined in equations (|19p - (|21|l can be written in terms of the new 
perturbation functions described by (|29|l and pO|) . as 

„-A 

e^^ai' {w cos 9 + w.e sin 6) + 2ai {v sin 6 — v^e cos 6)] , (39) 



Aur'^e sin ^ 



47rr-^e sin ^ 



« / 2 e' \ , 2ai 2A , . 2A . I -A • , 

<I> H 1 ai — e + 47re Ji > e it; sm f 

re 



I 2ai aie 
r e 



+2v ( ai' — — ^ — ^^^^ I cos S — e ^ai'w' sin S + 2aiii' cos i 



with <57i"'* = 0. 

Finally, the variations of pressure and density are described by 



5p 



5e 



IP 



--e w +vg) H — p 



/ -A / , \ , -A 

(— e w +v g) + . e w. 



(40) 

(41) 
(42) 



3.3 Boundary Conditions 

For axisymmetric polar perturbations, in general, the angular dependences of the functions w{t,r,8) and v{t,r,9) can be 
written as w oc sin (cos 6*) and v oc sin ^ 9eP^ (cos 0), where Pe{cos9) is the Legendre polynomial of order i. Thus the 
boundary conditions on the axis {9 — 0) will be 



ui = and v = 0, 

while on the equatorial plane (6* = n/2) we set 
w — and v^g = for odd £, 

w g — and v = for even £. 



(43) 

(44) 
(45) 



Finally, at the stellar center we impose the regularity conditions for w and v, i.e., w — and u = at r = 0. While the 
boundary condition at the stellar surface is that the Lagrangian perturbation of the pressure should be zero, i.e., AP = 0, 
which is equivalent to lu = and « = as we will use a barotropic equation of state and e = at the stellar surface for the 
background. 



4 NUMERICAL RESULTS 
4.1 Code Tests 

For our studies we adopted a polytropic equation of state (EoS) defined as 
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Table 1. The frequencies for /, pi and p2 modes, for the stellar model with M = 1.4Mq, for the harmonic indices £ = 2, 3 and 4 
in kHz. Here the superscripts (E) and (T) denote the frequencies found by solving the eigenvalue problem and those by evolving the 
2-dimensional perturbation equations, respectively. 



e 


/(E) 


/(T) 


(E) 


(T) 
p\ 




(T) 


2 


2.68 


2.75 


6.71 


6.87 


9.97 


10.21 


3 


3.26 


3.41 


7.65 


7.89 


11.07 


11.38 


4 


3.75 


3.89 


8.49 


8.73 


12.08 


12.38 



p = Knomb f — ) 
Vno / 



and e — nmt H , (46) 

no/ r - 1 

where nit and no are the baryon mass and normalized number density given by mi, = 1.66 x 10^^^ g and no — 0.1 fm~'^, 
respectively. Specifically, we adopt the values of F and K as F = 2.46 and K = 0.00936, whose values are in good agreement 
with the tabulated data for the realistic EoS A l|Pandharipandelll97ll ). For this EoS we focus only on the stellar model with 
mass M — 1.4M0, radius R — 10.35 km, and compactness M/R = 0.2. 

For the time evolution of the perturbation equations (|3ip , we use the iterated Crank-Nicholson scheme ( Teukolskv 2000l ) , 
which keeps second-order accuracy for the both of space and time. In order to eliminate the spurio us higher frequencies, we 
also add a 4th-order Kreiss-Oliger dissipation into the evolution equations IjGustafsson et al.lll995l ), where the coefficient of 
this dissipation term , ep, depends on the num e rical grid size. This type of numerical viscosity has been used successfully in 
Sotani et al.l (|2008al ) and lGaertig fc Kokkotas. I l|2008l ). Most of the numerical calculations in this paper were done in the 2D 
region of {r,9) with 50 x 40 grid points, where 6 is in the range of ^ 6 ^ 7r/2. The values of 9 = and 7r/2 correspond to 
the axis and the equatorial plane, respectively. It should be mentioned that the results were not significantly improved when 
a larger number of grid points has been used, i.e., 50 x 80 or 100 x 40. The numerical evolutions were stable for long enough 
to allow us to extract the specific frequencies for each spatial point by performing a Fast Fourier Transformation (FFT) of 
the time varying perturbation functions. Since in this paper we evolve the perturbed models for about T = 200 msec, in the 
FFT we can expect that possible maximum error in the estimation of the frequency A/, is around A/ ~ 1/T ~ 5 Hz. 

Before studying in detail the magnetar oscillations we performed test of the 2D evolution code on non-magnetized stellar 
models. In general for a spherically symmetric star which is not endowed with a magnetic field, the oscillations can be 
decomposed into spherical harmonics Yim which are independent of the angular index £. Thus by using the appropriate 
boundary conditions we can determine the specific frequencies by transforming it into an eigenvalue problem. The idea here is 
to compare the frequencies determined by the eigenvalue methods, to those derived from our evolution co de in order to check 
the ac curacy of the new method. For the solution of the eigenvalue problem we use the same code as in Sotani &: Kokkotaj 
1 20041 ). The regularity condition at stellar center, demands that the behavior of W and V in the vicinity of stellar center is 



of the form W ~ Cr^"*"^ and V ~ —Cr^/£ + ■ ■ ■, where C is some constant. Thus we set the initial data for the time 



evolution in the form 



w{r,e) 
v{r,9) 



er sin SPf (cos S), 
sin 9 dePi{cos t 



(47) 
(48) 



On a non-rotating cold neutron star without magnetic fields there exists only one family of fluid mode, i.e., the p (pressure) 
modes. The / mode is a fundamental mode of the pressure mode family, i.e., the po mode. There exists only one such mode, 
while there is an infinite number of higher p modes for every £. Traditionally the p modes are named as pi, p2, ... in ascending 
order of those frequencies. Note that the subscript used in the notation of the p modes shows in addition the number of radial 
nodes in the eigenfunction, while the eigenfunction of the / mode has no nodes. In table [1] we can compare the frequencies 
determined by the solution of the eigenvalue problem and those derived by evolving the 2D perturbation equations. We list 
results only for the /, pi and p2 modes with harmonic index £ = 2, 3 and 4. Since the oscillations with lower frequencies (or 
with longer wave-length) are not so sensitive to the number of grid points, the lower frequencies is expected to be calculated 
with higher accuracy. Actually, we find that the / mode with £ = 2 shows a deviation from the values derived by the eigenvalue 
problem of the order of 2.6% while the / mode with £ — 4 shows a deviation of about 3.7%. In any case the results derived 
with the time evolution code show a maximum deviation which is always bellow 5% even for the higher frequency modes. We 
consider that this accuracy is acceptable for the type of study that we perform here, i.e., to quantify the effect of the magnetic 
field on the fluid modes and to find out whether the polar Alfven oscillations are discrete or continuous. 



4.2 Polar Oscillations of Magnetars 

As initial data for the evolution of the perturbation equations we used the analytic forms for the perturbation functions w{r, 9) 
and v{r,9) described by the equations (|47p and (|48|) . We set initial data for angular indices £ = 2, 3 and 4 and we imposed 
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1000 2000 3000 1000 2000 3000 1000 2000 3000 



/ [Hz] / [Hz] / [Hz] 

Figure 1. The FFT of the perturbation function v{t,r,9) for a stellar model with B = 1 X 10^^ Gauss. The three lines in each figure 
correspond to three different angular positions, ~ 0, 7r/4, and it/2, and the different figures correspond to the different radial positions, 
r ~ 0, R/2, and R. In all figures we can see one peak for the Alfven oscillations, which is corresponding to 20-0 mode (i.e. the fundamental 
I = 2 polar Alfven mode), and another peak at higher frequencies corresponding to the frequency of the / mode with £ = 2 (2/). 




right panels correspond to the perturbation functions w{t,r,8) and v{t,r,9), respectively. 



at the equatorial plane the boundary conditions described by the equations (|44|l and (|45p . Actually, the perturbations on an 
axisymmetric background can be decomposed into two classes that satisfy different conditions at the equatorial plane, i.e., 
one class satisfies the condition specified in equation (|44p and the other one satisfies the condition given by equation H45p. 
These conditions describe perturbations for even and odd values of the angular index £. 

The various features of those oscillations are examined in two ways. The first involves checking of the FFT amplitude 
at various points inside the star. If the spectrum is continuous the peaks in FFT will depend on the position, while if the 
peaks are independent of the observer's position then they correspond to a discrete spectrum. The second check involves 
the study of the phase of oscillation for each peak frequency in FFT, i.e., for the discrete spectrum the phase should be 
constant throughout the star. Figure [T] shows the FFT of the perturbation function v{t, r, 6) at various points inside the 
star with _B = 1 x lO'^'' Gauss. The time-series has been generated by evolving an initial data set corresponding to ^ = 2 
deformations. In this figure the three different lines correspond to different angular positions i.e., S ~ 0, 7r/4, and 7r/2, while 
the different panels correspond to different radial positions i.e., r ~ 0, R/2, and R. It should be noted that the FFT of the 
other perturbation function, w{t,r,9) shows that the peaks correspond exactly to the same frequencies as in figure [1] and 
obviously the peaks are independent of the observer's position. The peak at 2754 Hz corresponds to the 2/ mode, while the 
lower peak at 300 Hz corresponds to the polar Alfven oscillation mode 20.0- Next we study the phase of the specific Alfven 
oscillation mode 20,0- In figure[5]we show the phase of the oscillation mode 2^0 for a magnetar with B = 1 x 10^^ Gauss. The 
left and right panels correspond to the phases for the perturbation functions w{t,r,6) and v(t,r,6), respectively. We observe 
that the phases for both functions w and v are almost constant except for some regions around ^ ~ 7r/3 for the function w 
and a part of r/R ~ 0.8 for the function v. We identify that these two regions with somehow strange behavior of the phase 
are related to the due to numerical loss of accuracy in the calculation of the effective amplitudes, where the eigenfunctions are 
almost zero because these are nodal points (see the left panels in figures |3] and |4| . Thus we can argue that the phase of the 
oscillation modes is independent from the position. With the above two tests we conclude that the Alfven oscillations of polar 
parity are described by a discrete spectrum, in contrast to the spectrum of the axial Alfven oscillations. This feature seems 
to be very similar to the case of inertial modes of rotating stars. That is, the spectrum of the axial inertial modes seems to 
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Figure 3. Effective amplitude of the eigenfunction for w{t, r, 8) for polar Alfven oscillation modes. The three panels correspond to Alfven 
oscillations with £ = 2 (left panel), with i = 3 (center panel), and with i = A (right panel). 




Figure 4. Effective amplitude of the eigenfunction for v{t, r, 8) with the several Alfven oscillation modes. The three panels are corre- 
sponding to the Alfven oscillations with £ = 2 (left panel), with £ = 3 (center panel), and with £ = 4 (right panel). 



be continuous, at least in the slow rotation approximation (e.g., Koiima 1 19981 ): Beyer fc Kokkotaa ( 19991 )). while the polar 
modes admit a discrete spectrum. 

Different initial data sets, such as initial deformation for £ = 3 and 4, verify the previous results, i.e., the spectrum is 
discrete and the oscillation modes have constant phase. In figures|3]and|31 we show the effective amplitudes of the eigenfunctions 
w{r,8) and v{r,9) with the corresponding peak frequencies in FFT, 300 Hz (for I — 2), 465 Hz (for £ = 3), and 500 Hz (for 
£ — 4), where we assumed that magnetic field strength is B = 1 x 10^® Gauss. 

In figure O we show the results of our examination for the dependence of the polar Alfven modes on magnetic field 
strength. This dependence of the polar Alfven oscillations on the magnetic field strength, can be eventually used to explain 
the observed higher QPO frequencies such as 150, 625 and 1840 in SGR 1806-20 or 155 Hz in SGR 1900+14. It should be 
noted that for weaker magnetic fields it is difficult to determine the frequencies of the Alfven modes because the corresponding 
peaks in FFT become a few orders of magnitude smaller compared to those of the fluid modes. This observation suggests that 
the polar Alfven oscillations will play no role on the dynamics of magnetars with weaker magnetic fields. In order to find the 
dependence of the polar Alfven mode frequencies to moderate and weak magnetic field strengths, one may study the problem 
by using mode analysis i.e. by d ecomposing the perturbations into spherical harmonics Yim- This type of study has already 
been done for Newtonian stars by Leel ( 2007 ). Actually, since we find that th e polar Alfv en modes form a discrete spectrum, one 
can safely study the problem by using modal decomposition. Additionally, iLed (|2007l ) observed that for oscillations confined 
in the crust there exists a critical value for the magnetic field strength. That is, the frequencies of the magnetic modes become 
smaller as the magnetic field becomes weaker and then, depending on the stellar model, they seem to disappear after a critical 
value. We have observed a similar disappearance of the modes, but we cannot give a conclusive answer since this might be 
a result of the numerical method that we used. Specifically, as the magnetic field becomes weaker the energy on the Alfven 
modes becomes minimal and the Fourier transform is dominated by the fluid modes, thus it becomes difficult if not impossible 
to extract the frequencies of the Alfven modes. Thus it is risky to conclude with time evolutions that there is a clear cut off 
frequency. Instead, in order to answer this, mainly academic, q uestion we plan to study the problem by using modal analysis 
of the perturbation equations . Concluding, the observation by Leel 1 200?! ) remains to be proved for the polar Alfven modes 
in the framework of general relativity. 

Finally, we examine how the fluid modes depend on the magnetic field strength. For a non-rotating magnetars with cold 
EoS there exists only the family of pressure modes (/ and pi modes) in addition to the polar Alfven modes mentioned earlier. 
However, these fiuid modes are hardly affected in any way from the presence of the ultra strong magnetic field unless its 
strength is larger than about 10^^ Gauss. This weak dependence on the magnetic field can be understood by realizing that 
restoring force for these fiuid modes is the pressure and the magnetic pressure is typically very weak compared to the pressure 
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Figure 5. The frequencies of the Alfven modes as functions of magnetic field strength. 
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Marks denote the corresponding frequencies determined by using FFT after time evolution, while the horizontal dashed lines show the 
frequencies of / and pi modes without magnetic fields. 





« 


:"4]?r"'""" 






• • 


15 


B [Gauss] 


= 2, 3 and 


4 as functions of magnetic 



of the fluid. In figure |6] we plot the frequencies of the / and pi modes with harmonic indexes £ = 2, 3 and 4 as functions 
of magnetic field strength. The marks correspond to frequencies of the modes of the magnetized stars while the horizontal 
dashed lines are the frequencies for the non-magnetized ones. From this figure, we can see that the effect of magnetic field on 
the frequencies of fiuid modes can be traced only when the magnetic field becomes stronger than a few times 10^® Gauss. Still 
the deviations between the frequencies for non-magnetized stars and those for magnetars with B — 4 x 10^® Gauss are very 
small such as 1.1 % for 2/, 1-3 % for 3/ and 0.3 % for 4/ modes while 0.4 % for 2P1, 0.3 % for 3P1 and 0.7 % for 4P1 modes. 
Since the frequencies of the polar Alfven modes depend on the magnetic field it is possible that for ultra strong magnetic 
fields to reach those of the / modes. 



5 CONCLUSION 

In this paper we study the polar type oscillations of strongly magnetized neutron stars. As a first step we derive the perturbation 
equations for an arbitrary magnetic field although for the numerical studies we adopted a stellar model with a global dipole 
magnetic field. In this study we ignored the effects due to the presence of the crust, since these types of modes are confined 
mainly in the core of the magnetar. By using the 2D time evolutions of the perturbation equations, together with the 
appropriate initial data sets, we estimate the Alfven oscillation modes as well as the fluid modes. We used two techniques 
to examine the form of the oscillation spectrum, i.e., whether is continuous as in the case of axial perturbations or discrete. 
Actually, by examining the FFT amplitude at various points inside the star and phase of the corresponding frequencies, we 
conclude that the Alfven oscillations with polar parity are described by a spectrum consisting only by discrete modes. There 
is a physical reason why the polar oscillations form a discrete spectrum while the axial ones a continuum. For the axial 
oscillations the restoring force is only the magnetic one and around the magnetic pole the direction of the magnetic force 
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becomes complicated, see also the discussion in ISotani et al] (|2008al ): IColaiuda et"all (|2009h : ICerda-Duran et all On 
the other hand, for polar oscillations the restoring force is not only the magnetic one but also the pressure. Thus in the absence 
of magnetic field these modes are still present while the axial ones reduce to a zero frequency spectrum. As we mentioned 
earlier, this behavior reminds the spectral properties observed the oscillations of rotating stars. There t h e axied modes that 
hadas restoring force only the Coriolis one were showing a continuous spectrum (see e.g. Koiima ( 1998l l: Bever fc Kokkota"^ 
while the polar ones were admitting a discrete spectrum . 

Finally the frequencies of the fundamental polar Alfven oscillations for typical magnetars are of the order of a few hundred 
Hz and their dependence on the magnetic field, can be used to explain some of the observed high QPO frequencies in SGRs. 

This is the first study of its kind for axisymmetric polar Alfven oscillations, it remains open to understand the form of 
the spectrum for non-axisymmetric perturbations while more co mplicated forms of th e magnetic field geometry should be 
assumed. For example, the inclusion of a toroidal component as in Sotani et al. ll2008b) w ill affect the form of the oscillation 
spectrum. Finally, the inclusion of the crust enriches the spectrum (jVavoulidis et al .l(|2008l)) and the various new features have 
to be taken into account in the attempts to understand the structure of magnetars via their QPOs. Actually, the introduction 
of crust, although it makes the problem more complicated might erase partially or completely the degeneracy in axial mode 
spectrum. 
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